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Abstract 


Recently,  Heading  and  Whipple  treated  the  problem  of  pro- 
pagation of  long  radio  waves  in  an  exponential  ionospheric  la- 
yer. They  dealt  with  the  case  of  an  obliquely  incident  plane 
wave,  a  vertical  magnetic  field, and  low  collision  frequency. 
The  case  treated  here  is  that  of  vertical  incidence,  oblique 
magnetic  field,  and  both  low  and  high  collision  frequency.  The 
results  obtained  do  not  agree  with  experimental  data.  This  is 
probably  due  to  the  fact  that  coupling  effects  were  neglected. 
An  exarrjple  is  given  in  the  Appendix  which  shows  that  the  omission 
of  coupling  terms,  hovrever  small  their  coefficients  may  be,  is 
not  always  justified. 
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_  1  _  '^ 

1.  Intro(Jiiction 

Recently  Hendin/'  rnd  tipple  trested  the  probler.  of  propagation  of 
lonp:  radio  wfives  in  an  ex^jonential  ionospheric  layer.  They  dealu  with 
the  cpse  of  pn   ©"blicjaely  incident  plane  wave  ^nd  a  vertical  niagnetic  field, 
and  er,suEed  both  v;ave  freqiiency  and  collision  frequency  to  be  smnll  cox- 
pared  to  gjTofrequency.  The  validity  of  the  latter  assumption  is  open  to 
question.  Jfievertheless  the  computed  results  are  reasonable.  The  overall 
reflectif^n  coefficient  is  found  to  increase  with  angle  of  incidence,  vhach 
agrees  with  exoerimental  data,. 

Because  o^  the  success  of  Heading  and  V/hipple  in  solving  the  problem  under 

the  ps«3un".ptions  mentioned  it  was  deemed  worthwhile  to  inveatipate  the  case 
of  vertical  incidence  and  an  oblique  magnetic  field.   In  this  esse  it  v;--s 
necessary  to  use  sonewhat  more  restrictive  assumptions  and  some  edditional 
approximations  in  order  to  calculate  the  reflections  coefficients.   Unfort- 
unately, th'^'  results  so  obtained  do  not  a^rree  well  with  experimental  date. 
In  viev:  of  this  fact  we  nay  wonder  Just  what  assumption  or  approyimation  is 
res-nonsible  for  the  discrepancy.   Since  the  approximations  used  in  this  paper 
and  in  that  o^f  HepdinP*  and  Whipple  are  very  similar,  the  similarity  between 
thp  results  of  both  papers  will  be  pointed  out  whenever  possible; 

2,  Oblioue  P-pfl^ction  In  Vertical  ?^>'g-netic  Field  Collision  Frequency  Sjall 

Ve  shall  reproduce  here,  in  bare  outline,  Heading"  and  Whipple's  theory 

of  obllnue  reflection  in  vertical  magnetic  field.   The  differential  e.;\ia.tinns 

2 

for   the   components   of   the   electric  fi^ld  for  this   case  were  derived  by  f<ilP:es    . 

They  are 

(2)  E  =  (m/iksine)dE  /dz 
where     m  = 


2n 
cos  ©  - 


1+iZ 


*  Kesults  were  computed  for  16  Tc  and  80  kc  waves  with  collislon-'l  frequency 
tJ  =  "5  T  10  sec.  ;  the  electron  density  was  assumed  to  hBve  Chapman  dis- 
tribution with  8  penetration  frequency  of  l.?7  mc  and  s   scale  height  of 
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2,  2 


03  -  wave  frequency, 

CO     -  critical   (plasma)   frequency 


Z 

k 


v/u 


coc 


2nA 


<j^-  -  gyro  frequency, 

V  -  collision  frequency, 
0  -  angle  betv/een  direction  of  incidence 
(in  the  xz-plane)  and  the  z  axis, 
and  the  time  dependence  was  taken  to  be  e    . 

If  Y  is  much  pireater  than  both  Z  and  unity,  and  X  increases  with  height, 
we  can  simrlify  the  problem  by  discussing  it  separately  in  two  regions. 
Region  I,  the  lower  one,  is  centered  at  the  height  given  by  X(z)  =  1.  In 
this  region  X  is  small  compared  to  Y,  thus  making  the  term  X/(l+iZ+YJ  small 
compared  to  unity.  In  Region  II,  the  upper  one,  centered  at  the  height 
given  by  X(z)  =  Y, 


z 

A 


t 


H 


Direction  of  incidence 
^x 


Oblique  incidence 
Vertical  magnetic  field 


z 


Direction  of 
incidence 


Vertical  incidence 
Oblique  magnetic  field 


X  is  much  greater  than  |l+iZ|,  and  tnerefore  m  is  small. 

In  Re-'ion  I,  then,  if  ^we  iron  the  term  X/(l+iZ+Y),  enuations  (l)  reduce 

to  two  separate  equations  for  E  and  E  .  Since  it  is  more  convenient  to 

X     y 

deal  with  E  than  with  E  we  can  make  use  of  equation  (2)  and  obtain  the 
z  X 

following  set  of  equations: 
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2 
^^^  -^  E  +  k^cos^eS  =  0 


^^^  ^(liiSE  )  .^E  =0. 

,  2   n   z    m   z 
dz 


The  latter  may  be  vritten  in  the  form 
^zz       1+m  "^     *    "^    m   z 


We  see  then  that  in  this  region  horizontally  and  vertically  polarized  waves 
are  propagated  independently.  The  horizontally  polarized  wave  is  unaffected 
by  the  medium,  as  is  evident  from  equation  (3).  If  we  substitute  the 
definition  of  m  into  (5),  we  obtain 

(6)  ^       ^   k^  (l^iZ)c°3^Q-X  g  =  0. 

^^      (1+iZ)  -  X 

We  have   here   a  wave  equation   for  ^,   \dx,h  the  refractive  index  changinp;  from 
cos©  to  unity  as  X  increases  from  zero  to  infinity.     For  small  values  of  Z 
the  refractive  index  has   typical  resonance  behavior  near  X  =»  1. 

If  X  increases  exponentially  with  height,   equatio:,    (6)    is  but  a   special 
case  of  a  wave  equation  treated  in  detail  ti^  Epstein     ,  if  we  Ipt   z^  -■  0   in 
,  the   e'lvu^tion  he   considereci    It  ■becomes 

V 

£  +  ^  e^ 

(7)  ^^^  *  ^     ^  '^'^y     5^  =  0   ,  r.  =  C(z-z^;  . 

1  +  e"" 

This  equation  can  be   solved  in  terms  of  hype rgeome trie  functions.      The  reflec- 
tion and   transmission  coefficients   can  be   obtained  in   terms  of  gar. ma   fimctions 
of  the  parameters  of  the  problem.     Whatever  the  variation  of  X  vri.th  hs:ight, 
the  vertically  polarized  component-  is  partially  reflected  at  thio  i-^vsl. 

In  Region  II  we  may  neglect-  t-he   second   term  in  equation   (1)    thus 
obtaining 

(8)  d^  .^.     ^     9       _  ^.,    X 


9 


dz' 


E     +  k  E     =  ik  ^  E 
X  X  Y     y 


Strictly  speaking  there  is  no  vertically  polarized  wavej  we  mean  here  the 
wave  polarized  in  the  plane  of  incidence. 
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(9)  d^  T7  ^  ,2   2^„     .,  X  ^ 

^  — TT  E  +  k  cos  Qu     =  -ik  V  ^  • 

These  a'"  -  coupled  wave  equations;  they  cannot  be  reduced  to  a  system  of 
independent  second-ordor  equations.  They  can  be  solved,  however,  in  terms  of 
generalized  hyperpeometric  functions,  if  X/Y  varies  exponentially  with  height. 
Four  independent  solutions  can  be  found,  T\<io   of  these  correspond  to  upgoing 
incident  waves,  and  the  other  two  to  dovmcoming  waves  incident  on  the  ionosphere 
from  the  top.  From  the  asymptotic  properties  of  the  former  tvro  solutions  we 
obtain  the  reflection  and  coupling  coefficients  for  Region  II, 

3,  Vertical  Propagation  in  Obliq-^e  Magnetic  Field;  Collision  Frequency  Small 

h 
R;/dbeck   has  shovTn  that  the  differential  equations  for  the  case  of 

vertical  inci'ience  and  oblique  magnetic  field  are 

,2 


(10,   fj  n,  M-  ^\  - 1  (ffi^)  ]  m 

dz         -        L        -  1+6         J 


i  TT       1-  (JklSl\       ^   --1  •  d5/dz        , 
"  "  2    I'-     dz   ^   -    ^2    ^   ■  ^     dz         -    .2 
+  1+0  x+o 

whPTP   thp   two  dependent  v^ripbles       |T^  and     \\  _  renlace   the   components 
E     Pnd  ■?!     of  the  electric  field  which  pre  used  directly  by  Wilkes.     Here 
the   totrl  field  is  first   resolved   into  componeiats   h?ivln,p:  polarization 
{V.   /E   )  eoTipl  u  pnd  l/n  : 

(11)  E;  =   (l-u^}-^(E^-uEy)    ,  eJ     =  uE;     , 


-1/r. 

y  "  X'    '  X  y 


E"  =    {1-v.^)      {t.   -uE    )    ,  E'     =  uE 

y  '         .y       X 


wnere 


(12)  u  =  i    [6  -   (62^1)^/2] 


(13) 


■    2^ 

1     sm  & 


2     cose         X  -  1  -  iZ 

2 
or,   in  terms  of  the  critical  collision  frequency  v^  =  w^sin  6/2cos0, 


Z 
(lU)  6 


"       X  -  1  -  iZ       ^ 
where  2^  =   '•?./">  . 


c        c 
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Then  we   define 


(15)  TT   =  (i-u'')^/^E, 


I  I     =   (1-u   )   '   t,„     . 

3.1  Coupling  'ferns 

We   see  that  the  coefficients  of  the   coupling  terms  in  equation  (10)    are 

(d5/dz)/(l+6   )   and  its  derivative,     V/e  have 

-Z 
f^y-s  d5/dz  c  dX 

1-.5"  z^  +   (X-l-iZ)^     ^^ 

c 

Ve  see  that  this  quantity  can  become  very  large  if  Z^Z  .  Instead  of  assuming 
that  V  «  61 T  ,  as  Whipple  does,  let  us  assume  here  that  v  «  v  .  This  coupling 
term  is  then  negligible  if 

(17)         dX/dz  «  kZ  . 

Since  Z  is  large,  this  condition  is  easily  satisfied.  In  the  remainder  of 
this  paper  coupling  terms  vill  be  neglected.  This  may  not  be  justified,  even' 
though  the  coupling  coefficients  are  snail.  It  may  be  that  the  lipiit  of 
solutions  of  the  equations  with  smaller  and  smaller  coupling  coefficients 
differs  from  the  solution  of  the  equations  with  coupling  terms  omitted.  An 
example  of  such  a  discrepancy  is  discussed  in  the  Appendix, 

3.2  Limiting  Polarization 

Since  we  are  interested  only  in  the  asymptotic  properties  of  wave  func- 
tions well  below  the  reflecting  regions  of  the  ionosphere,  we  need  to  consider 
only  the  ILmiting  polarization  (as  X  ->  0)  of  the  characteristic  waves  used 
here.  From  equation  (Ik)   it  is  evident  that  5  is  large  when  X  =  0,  Hence  u 
is  given  approximately  by 

-i     1+iZ 


(16) 


26     -2iZ 
c 


Since    |u|  «  1,  vre   are  dealing  here  with  two  approximately  linearly  polariz:?;i 
waves, 

3,3  Dielectric  Constants 

We   shall  now  investigate   the  dielectric  constants,   6_^,   in  equation   (10), 
These  are  given  by  the   Apple to n-Har tree  formula 


-   o   - 


(19)  e^     =1 2—1 "     u       U ~Vr 


2(X-l-iZ)   -  2 


L(X-l-iZ)^ 


+  Ul   cos  & 


In  the  IcA^rer  region   (Hericn  I  of  Heading  and  I\'hipple),  where  X  ^^  1     and 
therefore   X  «  I,   is^   obtain   the   foilovrinp;  approximations: 

(i)     For  the  ordinar7      (  +  )     ray  we  can  neglect  the  hY  cos  0     term  under 


the  radical.      The    (1+iZ)    term  is  negligible  compared  to 
Z ) .      Ke  nc« 

X(X-l-iZ) 


2        2 

Y  sin  &/(x-l-iZ).     Kence,  we  obtain 


(20)  £^-:^l  -     2       2         • 

Y  si  n't* 

(ii)      For   the   extraordinary   (-)   rpy  the   iiY  cos  Q  term  cannot  Tse 
neglected  'because  of   the   near  cancellation  of   the   radical 
and    the  preceding   term.     Ret?inin^   two  terms   in  the  'binoinial 
exppneion  of  the   mdic-l,  we  ^et 

(21)  z     -1-     ^^^""''^ 


2 

1+iZ-Xcos  fr 


Figure   1   shows  a  comparison  between  values  of  the   dielectric  constant 
obtained  from  £q.    (21)    (circles)    and   those    computed     from  the   Appleton- 
Kartree   formula  by  Kelso,  Thus  we   see   that  the    (+)    component  propagates 

alm.ost  as   if  the   ionosphere  were   not  there,   unless   sinQ  is  very  sm.all.      The 

other  component,   however,   behaves  q-.ite  differently!      the  dielectric  constant 

2 
shows  a   resonsnce,  and   it   changes   from  unity  at  X  =   0   to  sec  6     for  Ic  r5:e 

X.     "^he  dif ferr-ntial   ennrtion  satisfied  by   this   component   is 

2  2  2 

/^n\  d       TT       1  2        2,,      (l+iZ)cos  fr-Xcos  &     -n-         ^ 

(22)  — p-  Tl    +  k  sec  fr     -^- P I  I      =0. 

dz  ~  (1+iZ)   -  Xcos  & 

The   similarity  bet-^^en  this  equation   and  equation   (6)   is   quite   evident.      In 
the  case  of  exponential  layer 

x-x 

(23)  X  =  exp     -j^  J 

^f  Cfln  follow  F'pstein's  method  of  solving   the  equation. 

Let  2 

,^,  >                               -Xcos  &         ^^■'•iQ  • 
(^U)  u  =     =  e       ^^ 

1+iZ 


-  6a  - 


1.0'? 


F^ 
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:^h''ZOo 


0^-1 50  k 
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Figure  1 

CoiTiplex  Dielectric  Ccnstar.t  for  the 
Extraordinary  Kay  in  itepicn  I 

Approximate  Curi.'B 

O  E.'vact  values  ( ircn  Kelsc) 
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6  =  TT  -  taa'-'-S. 


Jiquation  (22}  becomes  then 

(25)      TT^      ^  u"^  n     +  kV    fu"^  +   (sec^O  -  lJ(l+u)-V^]  IT     •=  0    , 
(2^)  and  Epstein's   eq.    (l2)  nre   identical   if   in  the  latter  c     »  1, 


'1 
Co  =   f^»   P  =   0,      Hcj  we  put 

(26)  TT  =  u^^^f  , 
and  then  eq.    (25)  "becoaes 

(27)  u(u+l}f"   +   (2ikh+l)(u+l)r'   +   (ikh+ikh  sec6)(ikh-ikh  secO)f  =  0, 

or  the  hype rgeome trie  equation.     The  solutions  of  (25)    are   then 

(26a)  TT  "  =  u^  ^^  F(i.  ikh+ikh  secO,  i  ikh-ikh  secSj  t  2ikh+lj  -u), 

1 

or,   in  descending  powers  of  u, 

(28b)  TT  "  =  u'  ^^^  ^^^^  F(ikhiikh  sec6,   -ikhiikh  secOj-  2ikh  secft+lj-u""', 

2 

We  note  that  for  x  ->  oo ,   or    |u|   ->  oo  ,    the  TTo"  solution  is  appropriate, 
since  it  corresponds  to  an  outgoing  wave    (i-jith  the   time  dependence  exp(-iait), 
as  used  throughout  the  paper).     For  x  ->  -  oo  ,  or  u  ->  o,    the  TT^     ^'d  JJj^ 

functions  represent  incoming  and  outgoing  waves     respectively.     Their  limiting 
behavior  is 

(29)  TTi"  --  u"  '^^  =   (e^^'^)'  ""^  =  e*  ^^^^"^o^     e'^  ^^  . 

In  order  to  calculate  the  reflection  coefficient,  we  make  use  of  the  circuit 
relation 

(30)  TTo"  =  A(ikh-ikh  sece,-ikh-ikh  secft)   TT]_ 

+  A(-ikh-ikh  sec©,  ikh-ikh  sec©)  TTn" 

where 

r  „  N        r  (-<+p+i)  r(p--<) 
'^-•' '  r(-)  rw     ■  ■ 

We  have  here  ^n  erpression  relating  the  Bolution  of  the  protileni  ;,:h:  ch  sp.t- 
Isfies  the  outgoin.';  wave  condition  at  infinity  to  two  fimctions  which  ^eh^ve 
like  incident  and  reflected  wave  z.%    (-©).   OomTjinlng  (29)  and  (:^0),  we  obtain 


•       -  8  -  , 

thp  reflection  coefficient  as  the  following  ratio  of  the  complex  amplitudes 
of  the  incident  and  reflected  parts  of  the  total  wave  solution: 

(31)  vi  ^    n^i^h)  rC-ikh-ikh  sec&j    H-ikh-ikh  secQ  -^   1)   2khQ 
r\-2ikh}  P(ikh-ikh  secQ)  n(ikh-ikh  sec©  +  1) 


or 


-2 


-1. 


Pdkh-ikh  secQ)     "]      2kh(n-tan     Z) 


(32)       R  =     ^^^^^^      n^ikh) 

1-COS&     p(-2ikh)    [   n(-ikh-ikh  secft)  _ 

Using  the  relation 

2 

n±x)        =    pdx)    Pi-ix)   =  tt/s  sinh  vcc 

we   can  calculate   the  absolute  value  of  the   reflection  coefficient: 


(33) 


R 


sinh  Tikh(sec&  +  1)  L^    /      .     -1-7,1 
^ exp:2kh   (n-tan     Z) 

sxTih  nkh(sec6  -  1)  '-       •  -* 


^     1  -  e>:p[2nkh(l-secQ)J     ^^^^_,,^^  .^^h) , 
1  -  exp    [2nkh(-l-secOjJ 

Sj  nee   Z  is  a  function  of  freq;;ency,  define   Z     as  the   value   of  Z  when 


kh  =  1,  Then 
(3U) 


Z  =  v/ui   =  vch/kh  =  Z  /kh  . 


Equfitions  ('^'^)  nnd  i"^)  five   the  relation  between  the  frequency  and  the 
ms^Tiitude  of  the  reflection  coefficient  for  the  extraordinary  wave  in 
Region  I.  A  graph  of  Ie]  for  some  values  of  6  and  Z  is  shown  on  Fig.  2. 
In  the  upper  re,r:ion,.  nef^lecting  (l+iZ)  compared  to  X  or  Y,  we  get 


(35) 


e+  C::il  - 


^I^sin^G  +  Y  ,,    2^   Y^   .  ii^  xl/2 
—  ^ -F^  (Ucos  &  +  -^  sm  t>)  ' 

■  2r  ^^  r 


-1 


Evidently  the  dielectric  constants  are  functions  of  a  sinple  parameter 


(36) 


XA  =  .o^A..jj 


The  dielecLric  constarits  are  real   and  independent  of  collision  frequency,   just 
as  in  the  case   treated  by  VJhipple,     Equation  (35)   can  be    simplified   to 


N-> 


rn 

o 

o 
<; 

o 

o 

rn 

::;n 

o 

rn 
:^ 

rn 
o 


c^ 
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(37)      e^:s(l+  ^  tan^0)+  (r^sec^O+i  tan^Q)-'"/^ . 


From  this  we  see  that  e  ,  the  dielectric  constant  which  exhibits  a  resonance 
in  repiion  I,  is  positive  in  Rep:ion  II.  The  other  dielectric  constant,  e  , 
has  a  zero  at  r=l.  If  we  hpd  used  no  approximPtiona  this  zero  would  "he 
fo-'ind  ?t  r  =  1  +  Y~  ;  the  shift  to  r  =  1  is  due  to  our  having  neglected 
tprms  of  the  order  t~     comp^rpd  to  unity.  Figure  3  is  p  grpph  of 
e,  Pnd  c  for  9  =  "^0  .  ?or  lar^e  r  we  can  approximrte  c  "by 

+  mm 

12- 
(36)       e   =  1  +  -^  tan  8  +  rsecG  , 

while  for  small  values  of  r 

(39)      e^~l 

The   latter  values   are   identical  with  those   for  large   X  in  Region  I. 

In  Region  I  the  case  of  exponential  variation  of  X  wL th  height  can  be 
solved  exactly.      In  Region  II   the  problem  can  be   solved  approximately  by 
replacing   (37)   by 

(UO)  e_^  Si  (1+  I  tan^e)   +   (i  tan^©  +  rsecO). 

This  approxirription  hpe   thp  Pavrntp^e  of  being  ,e-ood  for  any  Q  for  "both  srr/ill 
pnri   Inrfp  values   of  r.     For  intermediate  values   of  r  it   is  good  for  sniall  -d. 
For  larrer  values   of  6  a   different  approximation  might  have   to  "be  used. 

In  Region  II,    then,  using  the   approximation   (UO)   and  exponential  varia- 
tion (23),    the  dielectric  constant  is  a  sum  of  a  constant  and   an  exponentially 
varying  term  for  both  waves.     VJe  can  consider  in  general  the  differential 
equation  » 

(Ui)  y_^  +  K^ey  =  0 


where 


e  =  1  +  expf  i«K   +   (x-x^)/h  J  . 


In  this  case  K   is  not  Lhe  free  space  wave  num.ber,  but  the  limiting  value 
of  the  v/Rve  mm'^er  as  x  — >  -  oo  .   If  wo  now  employ  the  ch.pnp-f   of  varir'hle 

(U2)         ^  =  exp  (x/2h)  , 

equation  (Ul)   transforms  into 


£ 
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£_ 


/ 


X 

0  0.5  rr 

e  IN     REGION    II 


I-'igure  3 
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(U3) 


1 


TT-     +  un  exD(ic^  -  - —  } 

d  *      n 


L  ^ 


y  =  0, 


This  eGuation  has  as  solutions  Bessel  functions  of  order  n  =  ^  2  i  t\  h  ^nd 

1        ^0 
argument  2hX  exp  ^  (i<K  -  r— ) •  Since  we  require  the  solution  to  behave  like 

an  outpoint;  wave  at  infinity, 


(liU)    y  =  H 


u( 


x-x 


"h  [2hexpl(i«.V') 


In  order  to  investigate  the  behavior  of  the  solution  at  (-oo)  we  express  it 
2iKj 

(1) 


in  terns  of  J„.  ,   and  J  o-  t. 
2iKh      -2rs:h 


(U$) 


n 


—. (  e     J  -  J  ) 

sm  nn         n    -n 


At  (-  oo)    the  argument  of  the  Bessel  functions  is  small  and  therefore 

P  (l+2iKh) 


„   ,^    r                   --.x-x      -i2iKh 
J^~e  h  exp   (   2"  +  -^ ) 


2h 


,  2iKh 
=  h  exp 


r  /  1    / 

p    [2Kh(n-'|)   +  iK(x-x^)j  /n(l+iKh)j 


J_^  --  [h  exp(|^ 


^      x-x       T-2iVih 

2h       ' 


P(l-2iKh) 


=  h"^^*^^  expr^<h»c-iK(x-x^)     /p(l-2iKh). 


The  reflection  coefficient  is  the  rai,io  of  the   complex  arrplitudes  cf 


exp 


r  *" 

and  exp     iK(x-x  ) 


Hence 


-iK(x-x^) 

(U7)      JRl    =  exp[-2)^h(n  -  1)   +  ■^h<'\   =  exp[-2<ch(n  -^)\. 

In  the  case  of  propagation  in  Region  II  we  have  the   special  cases  <=<^ 
and  o<     =  0 ,     From  ( UO )  we  have 

(U8) 
The  reflection  coefficients  are  then 


.2         ^        1   ^     2^  +  !■■        2_ 
i<         3  1  +  J  tan  ©  -  o"  tan  ft . 


(U9) 


"J   =  1 


exp(-2nKhsecS) 
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Thus  we   see   that  the  ~[\    wave   is  unaffected  in  Rerion  I  and  totally 
reflected  in  Rep-ion  II,  while   the    |  j     v:ave  is  partially  reflected  in  both 
regions.     For  reasonable  values  of  parameters,   the  partial  reflections  of 
the  77    wave  are   small.      Tnus  we  would  exnect  a  sin^jle  linearly  polarized 
component  in  the   dovmconinp  wave.      It  i?  well  known  that  at  very  low  fre- 
quencies a  sinfle   component  is  observed,     out  it  is  approximately  circularly 
polarized.      Hence  \<ie  have   to   question  the   validity  of  the   assumptions  made. 
The  most  vulnerable   is  the  one  concerninp:  the   smallness  of  collision  freq^aency. 
This  assumption  is  not  satisfied  in  those  regions  of  the   ionosphere  in  which 
the  reflection  of  low  frequency  waves  takes  place. 

U.     Oblique  Propagation  in  Vertical  Kagnetic  Field;     Collision  Freq..ency  Large 

V.'e   snail  now  investigate   the   case   of  large   collision  frequency,    Z  »,I. 

In  this  case,   in  Region  I  we  can  negl'-ct  X/  (l+iZ+Y)   in  equation  (1)    and, 

2 

in  -.ddition,  set  m  =  tan  &,  since  X  «  Z.  Thus  we  0Dta;n 

(50)  d^  „    ,  2   2„„    „ 
^  ^              — TT   iL  +  k  cos  Qh     =  0 

,  2  X  X 

dz 

(51)  -4^  ^k^cos^QE  =0, 

dz2  y  y     ' 

so  that  in  this  region  propagation  is  not  affected  by  the  presence  of  the 
ionosphere.  It  is  only  vhen  X  becomes  comparable  to  Z  that  we  get  an3'  effect 
at  all.  In  this  case  we  get 


(52)      f^[(^-)f^^x 


+  k^Cl+i  1)2^=0 


(53)      ^E^  k2(l.i|)  E  =0. 
dz   ^ 

In  the   limiting  case  of  normal  incidence   these   /^  auations   descrioe    the 

propagation  of  electromagnetic  vraves   in  a  medium  having    the   dielectric   constant 

6  : 

2 
(5U)  e  =  1  +  iX/Z  =  1  +  iu  /uiv   . 

Thus  we  find  that  in  this  approximation  the  medium  bscom.es  isotropic. 


^,     Vertical  Fi'cpapation  in  Oblique  J-iagnetic  Field'^     Collision  Frequency  Lcirtte 

Just  as   ir;   l.ht=   cose   of  low  collisioi:  frequency,    the   couplinj:  between    the 
v.-:o  cliarccteristic  waves    |  [      and  J]"     is  small  vjhen  v  »  ai^.     The  dielectric 
constants,    to   the   sarne  degree  of  approximation  as  used  in  Sec.   3.3, are 

(55)  e    ^1  +   DI/Z  =  1  +  ico  /ojv   , 

or  identical  with   (5U).     -""ron    (12)   and    (13)    It  is   readily  seen   that   the 
charrctpris.tictic  polrrizrtions  p.tp   circylar,  with  opposite  senao   of  rotation. 

Attain  the  medium  proves  to  be   isotropic  in  this  approxim.ation.      Yet  we 
knew  from  e^rperimental   dc-'ta  that  only  one  of  the   tvjo  circularly  polarized 
waves   is  reflected.      In  oropr  to  ©"btpln  this   effect  from   the   theory   it   is 
necessa.ry  to  retain  more   terns  in  the  die3ectric  constant  expression.     We  have, 
from  (19), 

(^^^  ^-^  -  l+i^Ycost,     ' 

the  we]l-known  longitudinal  propapation  approximation.     If  X  varies  as 
exp(z-z   )/h,   the  reflection  cocificients  are 

(57)  R^     =  exp(-2kh  tan--  ~^^)    , 

v:here  the   arctanfrent  is  restricted  to   the  range  0   to  k.     For  the  case   Z  »  Y, 

i 

vre  have 

f^r,.              ,     -1           Z                n       liYccs© 
^5fi)  ^^^^       l;t^co39    -  -? Z— 

The  ratio  of  the  reflection  coefficients  is  then 

(59)                   ^+  ~         /    i^v,  YcosS^ 
^  ~exp(-Ukh  — - — ). 

Thus  w-e  see  that,  if  the  increase  in  electron  density  is  very  gradual,  i,e,, 
when  kh  is  larpe,  the  ratio  of  the  two  intensities  is  very  snail, 

6,  Conclusions 

As  was  pointed  out   in  Sec.   3,    for  reasonable   values   of  parameters 
vp  v;o-ald   expect  p    sinfle  downconin^   linearly  polarl7.ed  component.      S-oeci- 
fically,   for  the   sane   values  of  parameters  as  used  by  Heading  and  V.Tiipple, 
vTith  the  electron  density  varying  as  exp   (z/0.6),    the  ratio  of  the   two  com- 
ponent,s  is  less  than  0.2  at  16  kc  a:id  less  than  0.01  at  60  kc.     Since   at 
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these  frequencies  the  dowiicoming  wave  usually  is  circularly  polarized, 
the  theory  apparently  does  not  apply.  It  is  hard  to  believe  that  this 
disagreement  could  be  due  to  the  difference  between  the  actual  electron 
distribiition  and  the  assumed  exponential  one.  It  is  more  probable  that 
the  neglect  of  coupling  and  the  assumption  of  low  collision  frequency  are 
responsible.  Neglecting  coupling  is  particularly  dangerous.   The  exai;:ple 
in  the  Appendix  shows  that  necrlec'oinp  coupling  in  a  particularly  simple 
case  can  lead  to  completely  vn'ong  results  with  respect,  to  transmission  and 
reflection  coefficients.  Although  the  example  used  is  a  rp'.her  special  case, 
it  is  by  -no  means  unreasonable,  and  the  equations  treated  then  bear  great 
resemblance  to  Rydbeck's  equations  of  propagation  in  the  ionosphere.  The 
agreement  between  the  experimental  data  and  the  high  collision  frequency 
results  is  better,  but  even  in  this  case  higher-order  approximations  have 
to  be  used  to  explain  even  the  most  fundamental  characteristics  of  propaga- 
tion. It  seems  evident  that  in  order  to  explain  fully  the  phenomena  of  low 
frequency  propagation,  the  theory  m,ust  not  assume  the  collision  frei:;uency  to 
be  larpe  or  suaall.  nowever,  if  these  assiw.ptionf.  ."re  drc.  ped  the  complexity 

of  the  eqiiatione  makee  it  doubtful  whether  the  pro'blem  can  "be  handled  cinalytlc- 
©lly.  A  nnmericsl  method  seems  indicated. 
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Appendix 
Let  us  consider  the   system  of  coupled  wave  equations 
(60)  u~       +(k     -  00  )  u~  =  2icA3u       +   iw     u* 

XX  X  X  J 

where  w  is  a  function  of  x  but  k  is  a  constant.      This  system  of  equations 

is  identical  with  Rydbeck's  equations  of  propagation  in  the  ionosphere 

(of.   Eq.    (10)),   if  we  let  e^   ■»   e_  -   1  in  the   latter  and   identify  co  with 

2   -1/2 
-(l/2)(d5/dz)(l  +  6   )     '      .       In  this  simple  case  we  can  solve  the  equations 

exactly.     We  define  a  new  set  of  dependent  variables  v" 


(61)  (0=f.'°^    '''"^)    {'''-] 

\v~'     Hsin©       cos©/     Vu"''     , 


where 


x 

(62)  0  =    r        a)(4)d4      , 

-co 


■L 


and  it  is  assumed  that  the  integral  exists  for  all  values  of  x,   and  that 
it  tends  to  a  constant,  ©   ,   as  x  ->  oo  ,     Eq.    (60)   then  transforms  into 

+  p     + 

(63)  v"    +  k     v"    «  0. 

Let  us  take  a  particular  solution  of  this  system,  say 


(6li)  V  "  e 


v~  -  0 
Inverting  (61),  we  have 


(65) 


u  \  /cos©  -  isin©\  /  v  \ 
-/"\^-isin©    cos©y  \  v  /  , 


Thus  we  see  that  the  asymptotic  forms  of  this  particular  solution  are 

+     ikx  +      ^   ikx 

(66)   x-^-oo^   u-^e    •      X-*  +oo  ^  u  -»•  cos©  e 

u"  -^  0   i  u'  ->  -isin©  e^^  . 

Hence,  so  long  as  ©  is  not  an  integral  multiple  of  2tt,  the  coupling  terms 
cannot  be  neglected.  If  we  did  neglect  them,  then  the  asymptotic  form  at 
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(-  00 )  vrould  be   identical  with  the   asyinptotic  form  at  (-!•  co  ) .     Furt.hermore 

•we  can  nake  the   coupling  coefficients  arbitrarily  sir.allj    and  yet,   if  ve 

keep  6     constant,    the  coupling  terms  never  become  negligible.     For  example, 

let 

9  =  8     e'^  (e-^'^  *  e'^)"^ 
(67)  ° 

ui  =  2a  S  (e  +  e      ) 

o  ' 

In  this  case,    aithour^h  0(oo)   =  S     is  independent  of  e,   the  extrema  of  co 
and  CO     tend   to  zero  as   e  ->  0.     Eeace  ve  see  that  it  makes  a  difference  at 

X 

what  stfige  e   is  allowed   to  approach  zero:   if  we  solve  Eqs.    (60)  nnd  then 
let  c  "$>  0  we  ^et   one  result,   and   if  we   let  e  ~>  0  in  2o.3.    (60)  and  now 
Bolvp  then,  we  ^et  another. 

The  aViove  es^mple  shows   clearly  that  the  omiseion  of  couplin^g  terms, 
however   small   the  coefficients  may  be,    is   not  always   justified. 
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Enrata 

Page  n       Eq.    (3.2)  In  denominator  of  inte^^ranrt   chanj^e 

«•    cos(:i/n  -  cos1(P+^)/tI   " 
to  read 

"    cos(ti/x)  -  cos[(i5+p')/T]  " 

Page  13,     5th  line  from       In  expression  for  lh-.(©,^)|    change 

"Ef"       to  read     «hJ" 


00  <^  lyOni 


Pa3e  1$,     0th  line  from       Change  "   pages    (I6)   and  (17/' 

bottom  ^^  ^^^  "pages    (17-a)    and   (16)   respectively" 


Pa-^e  IS,     Itth  line  from       Omit  "for  these   angles" 
bottom 

Page  15,     Uth  line  from       Change  "h(n,  ~)^     to  read  r'h^(7i,  ^) '  ■ 
bottom 

Page  17,      Change  "    |e(©,   p)  1  "    to  read   "  le-(9,   p-)  | "     in  caption  and  in 
designation  of  vertical  ordinate. 
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